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Abstract. We present in this paper a general model to determine the quantum 
properties of the light generated by a synchronously pumped Optical Parametric 
Oscillator (SPOPO) operating below threshold. This model considers time and 
frequency on an equal footing, which allows us to find new quantum properties, related 
for example to the Carrier Envelope Offset phase (CEO), and to consider situations 
which are close to real experiments. We show that, in addition to multimode squeezing 
in the so-called "supermodes", the system exhibits quadrature entanglement between 
frequency combs of opposite CEO phases. We have also determined the quantum 
properties of the individual pulses, and their quantum correlations with the neighboring 
pulses. Finally we determine the Quantum-Cramer Rao limit for an ultra-short time 
delay measurement using a given number of pulses generated by the SPOPO. 



1. Introduction 

nonlinear media are efficient generators of non-classical states of light, such as 
squeezed states or quadrature entangled states, which are valuable resources for quantum 
metrology and quantum information processing in the continuous- variable (C.V.) regime 
[U |2]. These non-classical states of light are produced by parametric down conversion 
either in the single pass regime |3] or inside a resonant cavity [4J. In the latter case, bright 
light is generated above some pump threshold value: one has an Optical Parametric 
Oscillator (0P0). Due to the resonant cavity buildup, OPOs produce non-classical 
states of light with the moderate pump power delivered by c.w. lasers. The quantum 
properties of c.w. OPOs have been extensively investigated over the past years, both 
theoretically and experimentally. In the non-degenerate case, they produce below and 
above threshold EPR entangled states of light (5j El Ej. In the degenerate case they 
produce below threshold highly squeezed light in the resonant frequency and spatial 
modes which are resonant with the cavity [8j [9]. 

Recently there has been a rising interest for OPOs pumped by trains of ultrashort 
pulses, because of their promising properties: firstly, with such pumps, the efficiency 
of the nonlinear process can be increased further and the oscillation threshold lowered; 



SPOPO time/ frequency quantum analysis 



2 



secondly the fact that the system generates light which spans over thousands of equally 
spaced frequency modes ("frequency combs") opens the way to the efficient generation 
of highly multimode non-classical light [TU1 ITT| IT2] ; in addition such sources can be used 
for the quantum metrology of ultra-short time delays [13J. 

Synchronously pumped OPO (SPOPO) are OPOs for which the cavity round trip 
time is synchronized with the time interval between the pump pulses in order to favor 
the resonant cavity-buildup. Such SPOPOs have been essentially developed as efficient 
classical sources of tunable ultrashort pulses [T4"t IT5l [T6| [T71 [I~8] in the femtosecond 
regime. They have also been used to generate squeezed light in the picosecond regime 
[19], and more recently to produce multimode squeezed light in the femtosecond regime 

Previous theoretical works on the quantum properties of SPOPOs below threshold 
[T2 | I2T1 122] have shown that SPOPOs generate multimode squeezed frequency combs and 
that there exist quantum correlations between pulses in different time intervals. These 
investigations have been carried out either totally in the frequency domain (longitudinal 
modes), or totally in the time domain. In this paper, we use the quantum input-output 
approach and develop a formalism which allows us to treat the quantum properties of the 
SPOPO in the frequency and in the time domain on an equal footing. This enables us 
to find new interesting properties, related in particular to the Carrier Envelope Offset 
phase (CEO), and to present a comprehensive model, capable of taking into account 
more realistic conditions, such as the intra-cavity dispersion [23] . 

The paper is organized as follows. In Sect. 2, we first present an input-output 
formalism for treating OPOs in the general case. Then, in Sect. 3, we treat the case of 
an OPO pumped by a train of synchronous pulses, and derive the quantum properties 
of the field that is generated below threshold. In Sect. 4, the quantum properties of 
individual pulses are derived. Finally we evaluate in Sect. 5 the quantum improvement 
of the estimation of ultra-short time delays using the pulses generated by the SPOPO. 
Appendices give more technical details on the derivation of some of the presented results. 

2. Input-output relation for lossless degenerate ring-cavity OPOs below 
threshold 

Let us consider an OPO in a ring-cavity configuration, where the cavity is transparent 
for the pump, and a x crystal in degenerate type-I configuration is placed in the 
cavity (see Fig. 1), Ml is the output coupler, and D represents the group of linear 
passive devices inserted for dispersion compensation. We will ignore the transverse 
spatial dependence for the sake of simplicity. In the free space outside the cavity, the 
positive-frequency components of the incoming and outgoing quantum fields can be 
written, within the slowly varying envelope approximation, as : 

E { +\t) = i<g a a (t)e- luJ(, \ with a = in,out, (1) 

where <r? is the single photon electric field, ojq is the carrier frequency of both fields, and 
ai n and a out are the slowly varying annihilation operators that verify the commutation 
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Figure 1. Sketch of the SPOPO ring cavity 



relations: 

[a a (t) , a a (t')]=0 and [a a (t) , a\ (t')]=5(t-t'). (2) 

One can also write the positive electric field operator in terms of annihilation operators 
of frequency modes, which are the Fourier transforms of time dependent annihilation 
operators and defined as: a(u) = J^a(t)e- iuit dt, so that: 

[a a (u),al(uj')]=2ir5(u-uj'). (3) 

Also for the sake of simplifying the model we assume that the losses are negligible. 
Thus, since the pump depletion is negligible below threshold, the 0P0 that we are 
dealing with is a linear lossless system, where the outgoing (output) field only explicitly 
depends on the incoming (input) field. Hence, the most general relation for our system 
reads, in the time as well as in the frequency domain: 



with ^11 (4) 






where C and S are integral transforms. For example, Ca in reads C(t, t')di n (t')dt' in 
the time domain and J^C_(u,u')a in (u!')du!' in the frequency domain, where C is the 
integral kernel [|J Since relation Q must preserve commutation relations ([2]), the above 
linear Bogoliubov transform should be complex symplectic (hereinafter, simply called 
symplectic), equivalent to a 2 x 2 real symplectic transform j25]. The necessary and 
sufficient condition for this is: 

n - l = (T l 'R}<T 1 , with O"! = ( J ^ V (5) 

Note that the effect of linear losses or other cavity configurations can also be described 
by using a symplectic transform, where the dimension of system is increased while the 
symplectic character of the transformation is preserved. 

For our linear and lossless system, the expression of 1Z can be found with the 
help of the relations between classical fields [8J. We denote Trt the cavity round-trip 

| In this paper, bold letters will denote matrix and vectors, and calligraphic letters will denote integral 
transforms. The underline will denote either the Fourier transforms, defined as: a(Lu) = a(t)e~ M *d£, 
or the kernels in the frequency domain. 
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transform, which relates the internal field envelopes after and before the output coupler 
Ml, i.e, 

where A\^{t) = E$ (t) exp(iuot) are the slowly varying envelopes. Since (6) also holds 
without the output coupler Ml, Trt is symplectic. We can write: 

Trt=Wib<SbaWai, (7) 

where «Sba is the transform of the crystal, and Mai and M±b are transforms from 
Ml to point A, and from B to Ml. These transforms must all be symplectic, and their 
explicit expressions will be discussed later. At Ml, the boundary condition for lossless 
coupling of the internal and external fields can be written as: 

(8) 

where t and r are real amplitude transmission and reflection coefficients (r 2 + t 2 = 1), 
which are assumed to be constant in the frequency range of interest. Then, with ^ 
and Q, we find: 

which, formally independent of the pumping condition, can be considered as a 
generalization of previous works (see, for example, Eq. (5.13b) of |2S])- From ([9]), 
it is easy to verify that 72-(7~rt) = ^- _1 ('7"rt)- Then, since Trt is symplectic, one can 
show with ^ that 1Z is symplectic. 

Outside the crystal, the transforms Mai and M%b, which are defined in ([7]), describe 
the lossless linear dispersive propagations. They are unitary diagonal, for example: 
Mai = Diag(WAi, W A1 ), where Uai is unitary and its kernel reads: 

U A1 (u,u') = 5(u - c/)exp[i0Ai(w)], (10) 

^AiO^-O being the spectral phase modification due to the propagation from Ml to A. 

In the lossless crystal, which extends from za = —l c /2 to z-q = l c /2, with the non- 
depleted pump approximation, the signal propagation can be described by the linear 
evolution of its slowly varying envelope: (A s (z), A* s (z)) T = S(z)(A s (za), A*(za)) t ■ So, 
the transform describing the propagation in the x^ 2 ^— crystal is given by: <Sra = <->(zb). 
From the classical description of traveling-wave degenerate parametric down conversion 
(PDC) pTH E7], S(z) can be found by solving 

d z S(z) = K(z)S(z), S(z A ) = 1, (11) 

with the kernel of K(z) given by: 
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In this relation the diagonal part describes the free propagation in the lossless medium 
and k s (—id t ) is the signal wavevector k s (ojQ + u) around Uq in which u has been 
replaced by —idt, whereas the anti-diagonal part describes the PDC process. Because 



criK\z)cri = — K(z), we have (TiS\z)cfiS(z) = 1. Thus, from Q, wee see that (11) 
and (12) describe a symplectic evolution. 

The effective nonlinear coefficient xo is given by: 

Xo = ( 3^3° , I d eS , (13) 

where uq = ck s (uo) / 'uq is the mean refractive index for the signal wave, srf C R is the 
effective area of degenerate three-wave mixing, and d e s is the effective nonlinearity 
coefficient for the given phase-matching configuration. The pump field envelope reads 
in the frequency domain: A p (z,u) = A p (u)exp[— ik p (u)z], where k p (u) is the pump 
wave-number at optical frequency 2uq + u. It is normalized so that |A p (2;,t)| 2 is the 
pump power. 

For OPOs, as long as the cavity finesse is relatively high, the amount of round-trip 
PDC is generally small below threshold. So, an approximate symplectic solution of (11 ) 
and (12) can be adopted, which is in fact the first order approximation of a Magnus 
expansion, also known as exponential perturbation theory, that preserves the symplectic 
character of the transformation 128 1: 



s B A=u BA g , (14) 

where Uba is crystal's unitary transform without pumping, and 




Q = exp I c n ) . (IT)) 
with 

S (u,u') =xdc£p{u + d)$(u,(J), (16) 

where = smA<f)(u,uj')/A(f)(u,u') is the phase-matching function with 

A(f)(u,u') = l c [k s (u) + k s (u') — k p (u + u/)]/2. Since Sq is symmetric, Qq is symplectic 
and positive-definite. Thus, from ^ and (14), we can write: Trt = Wrt£?rt, 
where Wrt = Uib^bpM ai is the round-trip unitary transform without pumping, and 
C?rt = U A1 G $Ll a_i is the positive transform of round-trip parametric gain. In the 
following we will replace C/rt by Q , which amounts to performing a unitary transform 
(a CT — > UxiCLcj) simultaneously on the input and the output. So, ^ can be rewritten as: 

n = r ~ l L Tj—n-- 1 )- 

V 1 - rWR T C/o / 

The round-trip unitary transform Wrt describes the passive round-trip spectral 
phase modification, 0rt(w). Around the signal carrier frequency uq, it can be 
approximated as a polynomial of u. For the sake of simplicity, we will assume the 
higher order dispersions are compensated or negligible so that 

0rtM = Ar T +wT o , (18) 
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where To is the round-trip time, and Art is the round-trip phase shift that detunes the 
resonance of cavity. Therefore, the kernel of U^r reads in the time domain: 

Urt(M') =6(t-t' + T )e i °i A ™, (19) 
which describes a time translation with a phase shift on the field envelope. 

3. Frequency combs as a basis for analyzing SPOPOs quantum properties 

We will now apply the general formalism of Section 2 to the specific case of SPOPOs, in 
which the pump pulses are synchronized to the cavity, so that the pump field envelope 
can be written as: 

A p (t) = ^4^e 2ifcAo a p (t - kT ), (20) 

k 

where S p is the pump pulse energy, 2A is the pump CEO, and a p (t) is the pulse envelope 
function, non-zero on ,% = [0, To] and normalized, such that f^. \a p (t)\ 2 dt = 1. Since 
the pump field is quasi T — periodic: A p (t + kT ) = A p (t)e 27,kA ° , we have a quasi periodic 



system. From ( 11 ) and ( 12 ), we see that the kernels of «Sba, Go and So are quasi periodic 
translation invariant, e.g., 

G (t + T , t' + T ) = e ifTlA °G (t, t')e- ifTlA °. (21) 



With (20), we can rewrite (|16j) in the time domain as: 

S (t,t') = J2Sc(t-kT ,t'-kT )e 2lkAo , (22) 

k 

where the kernel of S c reads in the frequency domain: 

S c (u,u / )= X0 l^ p 1/2 a p (u + u'Muj,uj / ). (23) 
In what follows, we will assume 

S c (t,t') = 0, Vt,t'<£%, (24) 

the validity of which requires that the spectral width of the phase matching function 
$(a;,u/) is much larger than the repetition rate T _1 , and that the pump pulses are 
localized around the center of each period and have a temporal width much smaller 
than Tq. These conditions are generally satisfied in practice with SPOPOs ultrashort 



pulses. Physically, (24) implies that the signal photons down-converted from the pump 



photons in each period will not travel into other periods when leaving the crystal. 



Due to the quasi-periodicity property (21), Qq commutes with a unitary pseudo- 
time translation IAq whose kernel reads: 

Uo(M') =5(t-t' + T )e- i(TlAo 

= U RT (t,t') e -W A *T+A ) ? (25) 
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where the second line comes from (19). So, Qq and Uq have common eigenvector 



functions. The eigenvalue problem of these two transforms is discussed in the Appendix 



A. With condition (24), we have: 



Uof£\e) = e»f£\e), 

with g n > 0, 9 e [— ir, n] and 

fn { , ) V2\±f*(-e,t) 

In the above expression, f n (9,t) is the eigenfunction of SoS^ and reads: 

1 +oo 



±9n f (±)/ 



/n(M) 



27T 



(26) 
(27) 

(28) 
(29) 



k=— oo 



with ^ n (t) = 0, Vt ^ where ^ n (t) is the eigenfunction of «S C «SJ with eigenvalue g\. 
The set {ip n } is orthonormal and complete over 

We see that each function f n (8,t) represents a continuous train of pulses, with a 
CEO of 9 + A . In the frequency domain, as shown on Fig. 2, f n (9,u) is a frequency 



/ 



Mi 



\ Y / / 

(2tm+A )/T 



Figure 2. Spectrum of a frequency shifted frequency comb. 

comb shifted by 9T$ l from the central frequencies of PDC {2mx + Aq)/Tq. The envelope 
functions {ijj n } are the so-called supermodes j21] and have also been studied in the 
context of single-pass and single pump pulse [[TDJ [TT], i.e., T — > oo. One contribution of 
this work is the introduction of the parameter of frequency shift, 9, that is proper for 
describing frequency combs. In Appendix B, we show that {f n {9,t)} forms a complete 
orthonormal set over the real time axis. So, the slowly varying envelope operators can 
be decomposed on the basis of frequency combs: 



a(t) = J2 



a n (9)f n (9,t)d9, 



(30) 



where a n {9) = f*(9,t)a(t)dt is the frequency comb annihilation operator, whose 
commutation relations read: [a n {9), a m {9')] =0 and [a n (9), a\ n {9')] =S nm S(9 — 9'). 

Now, on the basis of frequency combs, with ( 25p6 ), we are ready to perform a 



Bloch-Messiah reduction [29] on the transform (17). This yields: 

^in,n\9) 



&out,n{9) 



1 - re ie T n 



aL „ — C 



(31) 
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with 



Xn = ^(Art+Ao) f coshg n sinh# n . ^ 

sinhgn coshgn. 



We can rewrite ( |3l| ) under the form: 

a««t,«(0) = C n (9)a in>n (9) + S n {9)a\ nn {-9) (33) 



with |C n (0)| 2 — (^(O)! 2 = 1. Thus the input-output relation (33) describes a squeezing 
transformation if 9 = and a twin photon generation when 6^0. We conclude 
that SPOPOs are multimode squeezers for frequency combs with zero frequency shift 
from the central frequencies of PDC In the time domain, these frequency combs are 
continuous trains of pulses with a CEO half of that of the pump. For 9^0, SPOPOs 
generate pairs of entangled frequency combs with opposite frequency shifts, or in the 
time domain continuous trains of pulses whose CEOs are symmetric with respect to 
A , i.e., A ± 9. This is just as in the c.w. case where the c.w. modes with optical 
frequencies uq ± u> are entangled. These entangled pairs are in the so-called two-mode 
squeezed state, and each part alone of an entangled pair is in a thermal state [30J. 
This implies that, when detecting the squeezing of SPOPOs in a homodyne detection 
scheme, the local oscillator should be either a frequency comb with zero frequency shift 
or a combination of frequency comb pairs with opposite frequency shifts. 

From (30), we see that the resonant cavity enhances the squeezing effect in a 
frequency range having a width close to the cavity bandwidth. In practice, to obtain 
maximum squeezing the round-trip phase shift Art of the cavity should be adjusted so 
that Art + Ao approaches a multiple of tt. As can be seen from (31), the threshold 
of SPOPOs corresponds to the singular point of the input-output transform. Let us 
assume r > 0. Then, we find that the threshold of SPOPO is reached as 

1 + r 2 

g = max(#„) -> a cosh- — — — (34) 

2r cos(A RT + A ) 

for 9 = if there exists an integer n such that |Art + Ao — 2mr\ < tt/2, or at # = 7r when 
| Art + A — (2n + l)7r| < tt/2. This corresponds to the two regimes of the degenerate 
OPO above threshold HE1. 



4. Quantum properties of individual pulses 

Now, let us consider individual pulses. We define the annihilation operator for the pulse 
in the nth supermode and on the kth period, as: 

a n ,k = / ipn(t - kT R )a(t)dt 

JSTk 
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whose non-zero commutation relation is: a 



m,k' 



&nm$kk' 



line of (35), we have used (|29j) and (30). Inversely, one has 
a(t) = ^2a n ,k^n{t - kT R ). 



To obtain the second 



(36) 



n,k 



The Hermitian quadrature operators are defined as usual by x n ^ 

and p n>k 



t 

n.k 



2, with [x n> j,p m ,k] 



ibnmSkk 1 - From (31), (32) and (35), and 



assuming Art + Ao = 2mr, we obtain another form of the input-output relation for 
SPOPOs: 



l Jn,k 



s-l ±sg„A±) 
1 e %,n,k-s' 



(37) 



where q^l = x n ^ and q^l — p n> k- If ^rt + A = (2n + l)ir, another expression can be 
obtained by replacing the reflection coefficient r in the above expression with — r. From 



(37), we see that each output pulse can be considered as the combination of the direct 



reflection of a first pulse and of the transmission of the successive intracavity pulses 
related to pulses that previously entered the cavity. 

In what follows, we will assume that the input is in a coherent state, including 
vacuum state. Thus, the output pulses are in a Gaussian state characterized by their 



(±) 



variance matrix |25]. Using (37), the non-zero matrix elements can be found to be 



V 



t±) 



n,(j,k) 



with 
2 



r 2 + 



1 



2_ — ^2gi2gi n 
e ±2 9n 



2 1 — r 2g±2g n 



re 



±9n\\j~k\ 



(3J 



where Sj^ = 1 — Sj^, the first term of the right hand side is the quadrature variance, 
and the second is the covariance. 

Two conclusions can be drawn from this variance matrix. First, each pulse is in 



a Gaussian state which is not a pure state: one can verify that Vn(lk)^n(kk) > as 
long as g n >0. At threshold, the squeezing of the p— quadrature of each pulse is limited, 



r(-) 



(") 



^ — 2r 2 / (1 + r 2 ), whereas V^ k ^ — > oo. Second, there are pulse-to-pulse quantum 
correlations and anti-correlations on, respectively, the x- and p— quadratures. As the 
time difference between the pulses increases, these correlations decrease respectively at 
a rate of re ±9n . It should be pointed out, however, that these quantum correlations are 
not strong enough to ensure bipartite entanglement between any two pulses. This can 



be verified on the Duan criterion with the variance matrix given by (38) |31| . 



Let us finally consider iV successive pulses that are in the nth. supermode. As 
discussed above, these pulses are in a mixed Gaussian state described by the variance 
matrices of their quadratures ~V„ (N). It will be useful to look at the minimum 
eigenvalue problem of the p— quadrature variance matrix: 

1 2 



(39) 
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where a^(N) is the minimum eigenvalue, and c^g the corresponding eigenvector. 
With the explicit expression (38), they can be found semi-analytically [32j. The 



normalized elements of read: 



(N) 



e%j{N-2k-i) 



(40) 



where < k < N - 1, JC^ is a normalization coefficient, and 8^ can be found by 
solving: 

cos0g?(j\r + i) 



cos^^-l) " 
The minimum variance is then: 



re~ gn , with < N9^ < it. 



1 



(N) 



f g W n,0 Q~9 



(41) 



(42) 



It is interesting to note that, at large iV where — >■ 0, it becomes the variance of the 
^-quadrature of the squeezed frequency comb a n (0), and strong squeezing is re-found as 
expected. 

5. Quantum Cramer-Rao bound for time estimation using SPOPO pulses 

SPOPOs deliver trains of pulses which can be used as a clock when the pump repetition 
rate is locked on a time standard, with which can be measured with a very high accuracy 
the exact time of occurrence of a given event. It is well-known that quantum noise 
squeezing or quantum correlations can be used to improve the measurements of physical 
quantities. The question we address in this section is the determination of the minimum 
uncertainty that can be reached in the estimation of the temporal positioning of the train 
of pulses using the quantum properties of the light produced by SPOPOs, independently 
of the exact quantity measured on it and on the precise protocol employed for extracting 
the estimator from the measured quantities : this minimum uncertainty is given by the 
Quantum Cramer Rao bound [33j that we will now calculate. 

Let us consider a probe-field which is time translated by a small quantity r: 
a(t)e iuJot -> a(t + T)e iuJo{t+T \ which can be considered as the result of a unitary transform 



U (r) = exp r 



+oo 



a\t)(iuo + d t )a(t)dt. 



(43) 



We want to estimate the time delay, r, using N successive SPOPO squeezed pulses on 



the time interval [0, NTq]. Due to (36), the reduced system can be shown to evolve 

(44) 



under a unitary transform. In the Schroedinger picture, this reads: 



P iN \r) 



e inNT p {N) (0)e- inNT , 
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where p( N \0) is the initial state of the iV pulses, described by the variance matrix (38), 
and 

N-1 

= / ]^"mn,0'l n ,k^n,k, (45) 

k=Q m,n 

with n mn = j % ip^(t)(u - id t )4) n {t)<\t. 

According to the quantum estimation theory, the estimation precision of the time 
delay is limited by the Quantum Cramer- Rao bound, equal to the inverse of the quantum 
Fisher information F that depends only on the quantum state of the system used in the 



measurement, given by (44) in the present case. In Appendix C, we show that, under 
the assumption that the probe field is intense and r is small, F can be approximately 
evaluated as: 

FN = \Y}*' r f^ { -\N)a' n . (46) 

n 

where A„ (7V) = \Vn"\N)]~ 1 , and the vector elements of a' n are defined by: 

n.fc = y^^n,m(Q-m,fc), (47) 



a 



which have been assumed to be real. In order to minimize the Quantum Cramer-Rao 
bound with a fixed mean photon number in the probe-field during the measurement 



time, the vectors a' n should be proportional to 5 n ,oCoo\ where is defined by (39) 



Then, with (47), we find the optimal probe photon field envelope, which reads: 

N-1 

a<">(f) = a { Q N) ]Tcff (fc)^(t - kT R ), (48) 

k=0 

where ip' (t) is the probe-pulse envelope defined as: 

(u Q -id t )ijj' (t) = i>o(t), (49) 
and ckq = ^ Nno(ujQ + Aw 2 ), with Nuq = J \a^ N \t)\ 2 dt is the total photon number 
of mean probe field and Aw 2 , its spectral spread [13J. Then, the optimal quantum 
Cramer- Rao bound reads: 

Ar 2 = o J ff R (A0] A 2 2 , = 2Ar s 2 QI >^(A0] 2 , (50) 
2Nn (u% + Aw 2 ) sqllo,ov n , v i 

where At|q L is the standard quantum limit using iV coherent pulses, as given in |13| . 

Figure 2 shows an example of the quantum improvement of time estimation, 
Atsql/At, as a function of pulse number iV and for different levels of the distance 
to oscillation threshold. We see that this quantity is always larger than 1, as expected, 
and that the improvement increases with the number of pulses used in the measurement 



until it reaches a limit for large N equal to v / 2o'oo' ) ^ na ^ i s defined by (42). We can 

then know the minimum number of pulses, and therefore the minimum measurement 
time, needed to reach the limit. In addition, on can verify that the optimal quantum 



Cramer- Rao bound (50) is achievable with the optimal probe field (48 ) and the balanced 
homodyne detection described in ref [13j, where the local oscillator is composed of N 
equal intense coherent pulses in the supermode ipoit). 
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Figure 3. Quantum improvement of estimation of a time delay, At/Atsql, as a 
function of pulse number N. SPOPO cavity finesse : T w 30. 



6. Conclusion 

In conclusion, we have presented in this paper a novel time-frequency model for 
determining the quantum properties of OPOs below threshold based on a general input- 
output formalism that we have applied to the SPOPO case. The impact of the pump 
CEO and of the cavity round-trip phase on these properties has been studied. We have 
seen that the frequency combs, or continuous trains of pulses, form a natural basis for 
analyzing SPOPOs. As in the pulsed single-pass configurations, the squeezing of SPOPO 
is multimode, and the number of effectively squeezed modes depends on the pump 
spectrum and the phase-matching condition. We have shown that SPOPOs generate 
pairs of entangled frequency combs with frequency shifts symmetric with respect to the 
central frequencies of PDC. 

SPOPOs have also been analyzed in the regime of individual pulses. We have shown 
that each individual pulse is in a squeezed thermal state, and the two quadratures of the 
pulses are correlated and anti-correlated. The variance matrices of the quadratures have 
been analytically evaluated. Not surprisingly, due to these pulse-to-pulse correlations, 
strong squeezing is re-found in the multi-pulse regime. Finally, we have discussed the 
optimization of quantum improvement of time estimation with a limited number of 
SPOPO squeezed pulses. 

In the present paper we have applied the model to the simplest case of a lossless 
degenerate type-I ring-cavity in the dispersion compensated configuration, but it can 
serve as a framework for analyzing more complex configurations which are closer to 
real experimental conditions, such as higher-order intra-cavity dispersion effects or the 
influence of the transverse spatial dependence. 
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Appendix A 

Here, we first show how to find the common eigenvector functions of Uo = Diag(Wo, IAq) 
defined in (25) and 

Ho 



lnC?c 



So 

<Sn* 



(51) 



that is defined by (15). Since Uo is unitary, its eigenvalues are complex and on the unit 
circle, i.e., 

U f(9) = e ie f{6). (52) 
Writing f(9) = (fi(9)J 2 *(9)) T , from g), we have 
f 1 (9,t + T ) = f 1 (9,t)e^ Ao+d \ 

f 2 (e,t + To) = f 2 (6,t)e^- 9 \ (53) 

which means they are quasi T — periodic functions. Let us assume f(9) is also the 
eigenvector of Ho, i-e., Hof{9) = g{9)f{9), where the eigenvalue g{9) is a real number 
because Ho is Hermitian. From cr 1 'H cr 1 = —Ho, we see that a if (9) is also the 
eigenvector of Hq with eigenvalue —g{9). Moreover, it is easy to see that a if {9) is also 
the eigenvector of Uo- Thus, both f^ + \6) = f(6) and f^'{9) = a\f{9) are common 
eigenvectors of IA and Ho'- 

H f {±) (9)= ±g(9)fV(9), (54) 

Uof {±) {9) =e w f ± \9). (55) 
Now, since Hlf (±) (9)=g 2 (9)f {±) (9), with (53), we find: 

S S*f 1/2 (9)=g 2 (9)f 1/2 (9). (56) 



From (22) and (24), this means 



[ K 2 {t,t')f 1/2 (9,t')dt' = g 2 f 1/2 (9,t), 
J % 



(57) 



with K> 2 = SqSq and t e STq, where we have explicitly eliminated the dependence of 
g. Since IC 2 is Hermitian L 2 |3l], its eigenfunctions, {ip n {t)} defined by: 

° K 2 (t,t')i[> n (t'W = g 2 n Mt), (58) 



form a complete orthonormal set. Thus, with (55), (58) and (24), we find (29). 
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Appendix B 

It is easy to see that the functions of the set {f n (9,t)} are orthonormal. To show the 
completeness, one first notes that for an arbitrary well-behaved function, say a(t), its 
inverse Fourier transform can be written as: 

-+oo 

<\(t)-- — / «(w)e Mt du 

-oo 

i /•(fe+l/2)f2o+Ao/T 

— V / a(u)e iut du 
2tt k J(k-i/2)n +A /T 

-+n /2 

a B (u;,0e i(u,+Ao/To) Ma;, 

'-Q /2 

where f2 = 27rT ~ 1 , and 



i r 4 

1 

/ 

J-i 



a B (oo, t) = ^~Yl e im °Hkn + ou + A /T ) 



2tt 

k 

is the so-called Bloch amplitude, which is T — periodic. On the other hand, for each 
given 9, the set of functions {g n (9,t)}, defined as 

g n {9,t) = V^U9,t)e-^ e ^ t 

+oo 
k=~oo 

is T — periodic and is complete orthonormal basis for all T — periodic functions due to 
the completeness of {ip n }- Thus, we can project (Xb(w, t) on the basis {g n (uT , £)}: 

a B (u,t) = ^2a' n (uT )g n (ujT ,t), 

n 

with a' n (uT ) = f^, g^(uT 0: t)a B {oj, t)dt. Therefore, 

r+Oo/2 

a (t)= / V<( W T )^( W T ,t)e^ +Ao / T ^da; 
J-n /2 



To 



= /" *J]a n (0)/ n (0,t)d0, 



where we have made the substitution 9 — > uT in the second line, and 

a n {9) = ^—a> n {9)= / f:(0,t)a(t)dt. 

So, we have shown that any well-behaved function, for which a Bloch amplitude exists, 
can be written as the combination of the functions of {f n {&, t)}. Thus, frequency combs 
indeed form a complete orthonormal set. 
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Since our system, evolving as (44), is in a mixed state, the Quantum Fisher Information 
should be evaluated as |33j: 

F N = trpW(0)P N , (59) 

where Ln is the symmetric logarithmic derivative that verifies: 

^[L N pW(0)+pWL N ] = i[Cl N ,p N (0)}. (60) 

We can replace by A&n = &n ~ (&n) i n the above commutator. Assuming the field 
is intense and r is small, from (45), we have approximately: 

N-l 

AQ N = V2Y,Y1 a 'n,k A *n,k, (61) 
k=0 n 

where a! nk = J2 m ^n,m(S m ,fc) has been assumed to be real. 

Since the sub-systems of supermodes are uncorrelated, p( N \0) can be written as a 
direct product of sub-systems: p( N \0) = n^Lo ®Pn*\ where pn is the state of the nth 
supermode. It is convenient to project pn on the eigenstates of the p— quadratures 

m 

P { n N) =f \dW n (p,i?W\d N P<l N i?, (62) 

Jsfc2N 

with 

W n (p,j?)= [ d N xW n (x,^-)e~ iS ^\ (63) 

where W n (x,p) is the Wigner function. For our case, it is a Gaussian function of 2N 
variables. 

e -i(A£TAl +) (7V)A^+ApTA£-^(JV)Ap„) 

W n (x,p) = -T : j—. : , (64) 

I 2tt (AT) 1 1 1 2tt v£-> (AT) 1 1 

where An (iV) = [Vi (iV)] _1 , and Ax n = x — x „ and Ap n = p — p n , with x n and p n 
being the mean values. Now, since (p'\x\p) = —id p 5(p — p'), we have from (62): 

(p\[AZn,k,P ( n N) ]W) = i(d Pk + d P ' k )WMP)- (65) 
Then, it is straightforward to show that the symmetric logarithmic derivative reads: 

1 N ~ l 

^ = 7^EE <(U^K* A ^r (66) 

VZ n j,k=0 

Hence, the Quantum Fisher Information (57) now reads: 
1 N ~ 1 

i ^ = 2EE A £(ly)( JV X* a ^- (67) 

n j,k=0 
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